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The 2006 radio flare in the blazar CTA102 is the largest ever reported for this source. The analysis 
of the single-dish light curves revealed a possible shock-shock interaction as a scenario for the 
2006 outburst. In order to confirm this hypothesis we analyzed eight multi-frequency (2 GHz 
- 86 GHz) VLBI observations covering the 2006 flare. We detected a radio component ejected 
around 2005.9 and a stationary feature around 0.1 mas away from the core. We also performed 
a core-shift and spectral analysis on the fully calibrated VLBI maps. The frequency dependent 
position of the core slightly deviates from the typical v -dependence during the flaring state. 
This behaviour may reflect a mismatch between magnetic energy density and energy density of the 
relativistic particles. From the core-shift corrected maps we obtained the spectral parameters and 
computed the temporal and spatial evolution of the jet intrinsic parameters such as the magnetic 
field and the particle density. The overall picture of the source could be best interpreted as an 
over-pressured jet perturbed by traveling shock waves. 
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1. Introduction 

The technique of Very Long Baseline Interferometry (VLBI) offers the unique capability to 
resolve the structure of radio jets in Active Galactic Nuclei (AGN). Hence, we can study the mor- 
phology of those jets and trace the kinematical variations within long-term monitoring programs 
such as MOJAVE [Q]. Extending the single-frequency VLBI observations to multi-frequency ones 
provides additional access to the spectral evolution in AGN jets from which the magnetic field and 
particle densities are derived []5j The conditions in the direct vicinity of the black hole can 
be probed by analysing the frequency-dependent position of the innermost region of the jet, the 
so-called core-shift ||. 

The dense sampling of the 2006 radio flare in the active galaxy CTA 102 (z=1.037) with multi- 
frequency single-dish as well as VLBI observations provides an excellent laboratory to study the 
temporal and the spatial variation in the physical conditions of the source during the quiet and 
active states. The analysis and modelling of the single dish light curves revealed a travelling 
shock-recollimation shock interaction in a non-pressure matched jet as a possible explanation of 
the observed behaviour \M. Additional evidences supporting this hypothesis were found from the 
kinematic analysis of the long-term monitoring of the source at 43 GHz within the Boston Uni- 
versity Blazar Monitoring Program [||], at 15 GHz with the MOJAVE program [|J|, and from eight 
multi-frequency VLBI observations (2 GHz - 86 GHz) during the 2006 outburst. This study de- 
tected a new traveling component, which was ejected in 2005.9 with an apparent speed of (13±2) c. 
Furthermore, several stationary features along the jet were identified with the first one located at 
0. 1 mas from the radio core. The interaction between the travelling and the first stationary compo- 
nent occurred around 2006.3, which corresponds to the peak in the 37 GHz single-dish light curve 
and supports the hypothesis of a shock-shock interaction ph. 

Here, we present a core-shift and a spectral analysis on the fully calibrated multi-frequency 
VLBI observations of CTA 102 to further investigate the possible interaction between a travelling 
shock and a recollimation shock. Throughout the paper we define the spectral index, a, using the 
relation S v °< v a . The optically thin spectral index, Ofo, can be derived from the spectral slope, s 
of the relativistic electron distribution (N °< E~ s ), via the relation Ofo = — (s — l)/2. We define the 
optically thin spectral index as CCq < 0. 

2. Data and Data analysis 

We used eight multi-frequency VLBI observations from May 2005 until April 2007 span- 
ning a frequency range from 2 GHz to 86 GHz for our analysis. We used AIPS for the data cal- 
ibration and DIFMAP for the imaging and model fitting (see [||] for further details). Figure |] 
shows the uniform weighted VLBA clean images with fitted circular Gaussian components for 
the June 2006 observations at different frequencies. Since the absolute position of the source is 
lost during the data calibration, we used a 2D cross-correlation technique to reconstruct the rel- 
ative position among different frequency images taking into account the uneven (u,v) coverage 
between them. If we assume that the extended jet structure is optically thin, e.g., its location 
is frequency independent, we can use this structure to correct for the opacity shift of the inner- 
most regions |^ [KJ]. From the obtained shift and the position of the core at different frequencies 
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Figure 1: Uniform weighted VLBA CLEAN images with fitted circular Gaussian components at different 
frequencies for the July 2006 observation of CTA102. The lowest contour is plotted at 10 x the off-source 
rms at 43 GHz and increases in steps of 2. The observing frequency and the restoring beam size are plotted 
above each map. For the labeling we use capital letters for the same physical region in the jet and the numbers 
increase with inverse distance from the core. We adopt the following cosmological parameters: £l m = 0.27, 
£2a = 0.73 and Hq = 71kms~ 1 Mpc~ 1 . This results in a linear scale of 8.1 1 pcmas -1 or 26.45 ly mas ~' for 
CTA 102 (z=1.037). With these conventions, lmasyr -1 corresponds to 52.9c. 



we calculated the core-shift and derived the magnetic field and particle density following []7j, ^J. 
Once the images were corrected for the opacity shift, we convolved them with a common set of 
image parameters (i.e., beam and pixel size), stacked them together and fitted a homogeneous 
synchrotron self-absorbed spectrum to each pixel to obtain the turnover frequency, v,„, the turnover 
flux density, S m , and the optically thin spectral index, Oo- A homogenous synchrotron self-absorbed 
spectrum is given by S v « S m (v/v m ) 5 fl -exp ^-T m (v/v,,,)" 0-5 / 2 ^ / (1 -exp(-T m )), where 

T m ss 3/2 ^(1 — 4ao/15) 1 ^ 2 — 1^ is the optical depth at the turnover. The uncertainties of the ex- 
tracted turnover values were computed from Monte-Carlo simulations taking into account the un- 
certainties of the image alignment and the SNR-based flux density error in each pixel. 



3. Core-shift and spectral analysis 

The frequency-dependent core position can be approximated by r core = A ^v -1 ^' — v r J^ kr 
where v re f is the reference frequency, and k,- reflects the physical conditions in the source, i.e., 
the evolution of the particle density and the magnetic field along the jet 1 [|7|, g]. We applied 
this approximation to the measured core-shfits and obtained an average value of (k r ) = 0.8 ±0.2 
and (A) = (3.4 ± 1.6) mas/GHz. Using the Equations presented in [|7| |8]] and assuming k r ~ 1, 

L £ r = ((3 — 2ao)b + 2n — 2)/ (5 — 2ao) where b parametrizes the evolution of the magnetic field B °= R~ b and n the 
evolution of the particle density Af °= R~" . If there is equipartition between magnetic energy density and energy density 
of relativistic particle, k r = 1 for b = 1 and n = 2. 
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CXq = —0.5, and a lower electron Lorentz factor at core of y m i n , C ore = 100, we computed the av- 
erage magnetic field and particle density at the 22 GHz core to be (5 C ore,22GHz) = 100mG and 
(Acore,22GHz) = 60cm -3 (see ^ for more details). In Fig. || we present the temporal variation of 
^core,22GHz and A^ CO re,22GHz derived with the assumptions mentioned above. 
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Figure 2: Evolution of the core magnetic field, B cor e (top), and the relativistic particle density at the core, 
Ncoie (bottom), as derived from the core-shift analysis. 

For a distance r < 2 mas from the core, the spectral turnover was within our frequency range 
and we were able to extract v m , S m , and do. However, farther downstream, the turnover frequency 
moved out of our frequency range and we could only extract the spectral index, a (S v °< v a ). For 
r < 2 mas, the turnover frequency decreased from v m ~ 40 GHz to V ~ 5 GHz with a local maximum 
of v ~ 8 GHz at r = 1.5 mas. A similar behaviour was found for the turnover flux density, declining 
from S m ~ 2.2 Jy to S m ~0.2Jy with a second peak at r = 1.5 mas of S m = 0.6Jy. At larger distances 
from the core we found additional local maxima in the evolution of a at r ~ 4mas, r ~ 9 mas, and 
r ~ 16 mas. In Figure || we show the temporal evolution of a along the jet for 2 mas < r < 10 mas. 

Based on theory of synchrotron self-absorption we used the detected spectral turnover for 
r < 2 mas and derived the magnetic field and particle density along the jet [^, ||]. The magnetic 
field decrease from B ~ 60mG at r = O.lmas to B ~ 2mG at r = 2.0mas. For r < 0.6mas the 
magnetic field geometry was best described by a toroidal field (B oc R~ - 9 ) and farther downstream 
by a poloidal one (B °c R~ 2 ). For the calculation of the particle density we assumed both adiabatic 
and radiative losses and computed the evolution of the upper and lower electron Lorentz factor 
along the jet using the derived evolution of the magnetic field. The particle density took values 
between 100cm -3 and 0.8 cm~ 3 with a secondary maximum of 5cm~ 3 at r = 1.5 mas. Based 
on the magnetic field and the particle density we derived the magnetization along the jet a mag = 
Ub/U e , with Ub = B 2 / (87r) the magnetic energy density and U e the energy density of the relativistic 
particles. The magnetization decreased from CJ mag ~ 0.1 at r = O.lmas to 0.001 at r = 

2.0mas. At r = 1.5 mas we computed a second peak in the magnetization of 0.03. 

4. Discussion and Conclusions 

The flaring activity in blazars is in general explained with the propagation of a relativistic 
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Figure 3: Evolution of the spectral index, a (S v <x v a ) for (2mas < r < lOmas) along the jet ridge line (see 
[ |I"2| l for details on the jet ridge line). For reasons of readability only one error bar per epoch is shown. Left: 
values along the jet ridge line. Right: average values transversal to the jet ridge line. 



shock wave through an underlying steady state jet. At the shock front, the magnetic field and 
particle density increase, which leads to an enhancement of the emission (see, e.g., [11]). Due to 
the limited resolution, we first detect this effect in the increased core flux density. From the analysis 
of the single-dish light curves we derived that the 2006 radio flare started around 2005.6 [|J], which 
corresponds to the first increase in B core and A^ore- In addition, we detected a traveling component, 
associated with the traveling shock around 2006.3 [0], which is in agreement with the decrease in 
Score and A^core- Thus, we interpret the variation in the magnetic field and the particle density to 
be due the propagation of a travelling shock through the core region. The second increase in B core 
and A/core can be explained with a new traveling shock entering the core region (see Fig. 17 in [@|). 
The increase in the spectral index, a, along the jet can be interpreted by i) traveling shocks or ii) 
recollimation shocks. If the increase is caused by a moving shock, the position of the peak in a 
propagates downstream and is stationary in the case of a recollimation shock. In the latter case the 
increase in a is caused by an adiabatic compression of the plasma at position of the recollimation 
shock [14]. Based on this differentiation and on the results of the kinematic analysis [Q], we assign 
the increase of the spectral index at r ~ 4 mas and r ~ 16 mas to recollimation shocks. Since the 
increase in a around r ~ 9 mas is slightly traveling downstream it could be caused by a traveling 
shock wave. The increase in the turnover values at r ~ 1.5 could be the result of an interaction 
between a traveling shock (j3 app ~ 8.5c) with a recollimation shock at r ~ 1.5 mas (see also Fig. 
8 in [Q]). The temporal increase in the spectral index for r < 0.5 mas indicates the passage of a 
travelling component through this region. However, the limited resolution does not allow further 
investigation within that region. 
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The derived magnetic field and particle densities from the spectral analysis are in good agree- 
ment with the ones calculated from the core-shift analysis. The obtained value for the magnetiza- 
tion of a mag < 1 reflects the slight departure from k r ^\ and points towards a particle-dominated 
jet. The change in the orientation of the magnetic field from toroidal to poloidal could reflect a 
spine-sheath jet, i.e., for r < 0.5 mas we observe mainly emission from the spine and for larger 
distances emission from outer sheath. 

Summarising, we derived the magnetic field and particle density of the jet in CTA 102 during 
the 2006 radio flare and the obtained variation is in agreement with the passage of a relativistic 
shock through the core region. Based on the combination of the spectral analysis presented here, 
and the results of the kinematic analysis in [Q] we claim that the observed behaviour in CTA 102 is 
the result of the propagation and interaction of traveling shock waves with recollimation shocks in 
an over-pressured jet. Further analysis of the interaction between these features will be conducted 



using relativistic hydrodynamic (RHD) and emission simulations following J jl% |14| ]. 
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